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A confined incompressible clastic film does not deform uniformly when subjected 
to adhesive interfacial stresses but with undulations which have a characteristic 
wavelength scaling linearly with the thickness of the film. In the classical peel ge- 
ometry, undulations appear along the contact line below a critical film thickness 
or below a critical curvature of the plate. Perturbation analysis of the stress equi- 
librium equations shows that for a critically confined film the total excess energy 
indeed attains a minima for a finite amplitude of the perturbations which grow with 
further increase in the confinement. 
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1. Introduction 

Spontaneous surface and interfacial instabilities of thin liquid films have been stud- 
ied in different contexts, e.g. the classical Saffman- Taylor (Saffman et al 1958) 
problem in a Hele Shaw cell in which flow driven fingering patterns develop at the 
moving interface of two viscous or viscoelastic liquids (Homsyet al 1981, Nittmann 
et al 1985); disjoining pressure induced rupturing and dewetting (Reiter et al 1992, 
Sharma et al 1998) of ultra thin viscous films; spiral instabilities in viscometric 
flow of a viscoelastic liquid (Muller et al 1989, McKinley et al 1995); and finger- 
ing instability and cavitation during peeling a layer of viscoelastic adhesive (Fields 
et al 1976, Urhama et al 1989). While most of these viscous and viscoelastic sys- 
tems have been well characterized experimentally and theoretically, similar surface 
undulations of confined thin elastic films pose a different kind of problem despite 
geometric commonalities with many of the liquid systems. The essential difference 
being that unlike in the liquid system there is no flow of mass and consequent 
permanent deformation in the elastic body, where the extent of deformation is gov- 
erned by the equilibrium of the external surface or body forces on the material and 
the elastic forces developed. 

The specific system that will be described in this paper is a thin layer of elas- 
tic adhesive confined between a rigid and a flexible plate. While the film remains 
strongly bonded to the rigid substrate, the flexible plate is detached from it in 
the classical peel geometry. High aspect ratio of such systems is noteworthy as the 
lateral length scale far exceeds the thickness of the film resulting in high degree of 
confinement for the adhesive. As a result, adhesive stresses at the interface does 
not always result in uniform deformation through out the whole area of contact, 
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rather spatially varying deformations (Gent et al 1958, Gent et al 1969, Ghatak 
et al 2000, Monch et al 2001) attain lower energy for the system. Experimentally 
we see the existence of a critical thickness of the film or a critical curvature of the 
flexible plate below which the contact line between the film and the flexible plate 
does not remain straight, but turns undulatory with a characteristic wavelength 
which increases linearly with the thickness of the film (Ghatak et al 2003). While 
experimentally this phenomenon has been characterized, there is not much under- 
standing as what drives this instability in a non-flow purely elastic system and how 
does the curvature of the plate or the thickness of the film result in critical con- 
finement of the film. Here I present a perturbation analysis which addresses these 
questions highlighting the duel effects of the incomprcssibility of the elastic film 
and its confinement. 



Figure 1. (a) Schematic of the experiment in which a model elastic adhesive remains 
bonded to a rigid substrate and a flexible plate is detached from it with the help of a 
spacer inserted at the opening of the crack. For a critically confined film, the contact line 
does not remain smooth but becomes undulatory as shown in video-micrograph (b). 

The schematic of our experiment is represented in Figure la in which an elastic 
film of thickness h and shear modulus /U remains strongly bonded to a rigid substrate 
while a flexible plate of rigidity D in contact with the film in the form of a curved 
elastica is supported at one end using a spacer of height A. The straight contact line 
between the film and the plate becomes wavy when the film thickness h decreases 
below a critical value h c or the curvature of the plate decreases below a critical value 
l/pc- Figure lb represents a typical video-micrograph of such undulations which are 
characterized by two different length scales: the separation distance A between the 
waves which scales as A ~ Ah and the amplitude A which varies with D and /i as: 

1 /3 

A <~ (D/fi) ' (Ghatak et al 2003). The figure depicts also the co-ordinate system 
in which x, y and z axes represent respectively the direction of propagation of the 
contact line, the direction of the wave vector and the thickness coordinate of the 
film. The y axis is located along the tips of the waves so that the film is completely 
out of contact with the plate at < x < a. Assuming the film to be incompressible 
and purely elastic with no viscous effect, we write the following stress equilibrium 
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relations in the absence of any body forces, 

Px {^xx ~t~ V"yy ~l~ l^zz) 

Py = V (Vxx + Vyy + V zz ) 

Pz = M {Wxx + W yy + W zz ) (2.1) 

where, u, v and w are the displacements in the x, y and z directions respectively 
and p is the pressure. Here and everywhere s x = ds/dx and s xx = d 2 s/dx 2 . The 
incomprcssibility of the film results in: 

u x + v y + w z = (2-2) 

These equations are solved using the following set of boundary conditions (b.c): 

(a) since the film remains strongly bonded to the substrate we use no slip boundary 
condition at the interface of the film and the substrate (z = 0), 

u(z = 0) = v(z = 0)=w(z = 0) = (2.3) 

(b) we assume frictionless contact at the interface of the film and the cover plate, 
z = h which results in zero shear stress at the interface, 

o- X z(x,y,h) =0 = a yz (x,y,h) (2.4) 

(c) we assume continuity of normal stress across the interface (z = h) which implies 
that the pressure is equal to the bending stress on the plate 

p(z = h) = DV 2 iP(z = h) at s<0 (2.5) 

here V = d 2 /dx 2 + d 2 /dy 2 is the two dimensional Laplacian, D is the flexural 
rigidity of the plate and tp = w\ x z=h is its vertical displacement. Since the plate 
bends only in the direction of x co-ordinate its vertical displacement tp remains 
uniform along the y axis; hence, we simplify b.c. 2.5 as 

p (z = h) = Dtp xxxx (z = h) at x < (2.6) 

At < x < a there is no traction either on the film or the plate which yields 

Vxz\ z=h = Gyz\ z=h = Ozz\ z=h =0 (2.7) 

Equation (2.1)-2.2 can be written in dimensionless form using the following dimen- 
sionless quantities, 

X = xq, Y = y/h, Z = z/h, U = uq, V = v/h, W = w/h, * = tp/h, P = pe 2 /^ 

While thickness h of the film is the characteristic lengths along the y and z axes, q^ 1 
is that along x. The length q^ 1 can be derived as the ratio of the deformability of the 
plate and the film (Dillard et al 1989, Ghatak et al 2004): q- 1 = (Dh 3 /3/j) 1/6 . The 
quantity e — hq defined as the ratio of the two characteristic lengths is a measure of 
the confinement of the film such that a lower value of e represents a more confined 
film. Equations 2.1 and 2.2 can then be written in the following dimensionless form 

Px = t 2 Uxx — Uyy + Uzz 

Py = e 4 Vxx+e 2 (V Y Y + Vzz) 

P z = e A Wxx+e 2 (W YY + W Z z) 

= Ux + Vy + Wz (2.8) 
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while the boundary conditions 2.3-2.7 results 

(a) U(Z = 0) = V(Z = 0) = 

(b) a xz (X,Y,Z=l) = = <T YZ (X,Y,Z=l) 

(c) P(Z = 1) = W XXXX atKO 

(d) = ^ XXXX at < X < aq (2.9) 

where aq is the dimensionless crack length. Equation 2.8 is solved by the regular 
perturbation technique which assumes that the solutions consist of two components: 
the base solutions which remain uniform along the Y co-ordinate and the correction 
term which incorporates the spatial variation along the Y axis. Thus the base 
solutions are of order e° and the perturbed solutions are of order e 2 , e 4 , . . ., so that 
any variable T {X, Y, Z) = T a (X, Z) + e 2 Ti (X, Y, Z) + e 4 T 2 (X, Y, Z) + . . . where 
T = P, U, V and W. Inserting these definitions in equation 2.8 and separating the 
base (Y independent) and the perturbed (Y dependent) terms yield: 

P 0X = {e 2 U oxx + Uozz) , Poz = (t 4 W oxx + e 2 W 0ZZ ) , = U ox + W 0Z (2.10) 
and 

e 2 P lx +e 4 P 2X = e 2 (U 1YY + U lzz )+e i (U lxx + U 2YY + U 2ZZ ) + e 6 U 2XX 
t ny -t- e 4 P 2Y = e 4 (Vi YY + Vizz) + e 6 {Vi xx + V 2YY + Vizz) 
e 2 P lz + e 4 P 2z = e i (W 1YY + W lzz ) + e (i (W lxx + W 2YY + W 2zz ) 

= U lx + V 1Y + W lz (2.11) 

which are solved using b.c. derived from equations 2.9a-d. 

Base Solution: Since for a thin film, e 2 << 1, equation 2.10 can be simplified by 
neglecting the terms containing e 2 . Integration of the resulting equations (presented 
in detail in reference Ghatak et al 2004) finally leads to the following solution for 
the base components of the displacements in the film and the plate, 

U = (3Z 2 /2 - 3Z) F'fa (X) , Wo - (3Z 2 /2 - Z 3 /2) F'fo (X) 

0i = e x/2 (aqe x/2 + (3aq + 4) sin (VSX/2^ /V3 — aq cos (V3X/2)) 

<j) 2 = e x/2 (aqe x/2 + (3aq + 2) sin (V3X/2) /V3 + (aq + 2) cos ) 

F' = 3A/ (6 + I2aq + 9 (aq) 2 + 2 (aqf^j 

*o = F'(f> 2 (X) X < 

= F' (2 (aq + 1) + (3aq + 2) X + aqX 2 - X 3 /3) < X < aq (2.12) 

which suggest oscillatory variation along X with exponentially vanishing amplitude 
away from the contact line. Such displacements results in the following expression 
for the dimensionless work of adhesion (Ghatak et al 2004) G = Wa/ (p-/q)'- 



G = g (aq) (27 A 2 /2e(aqy 
g (aq) = 8 (aqf (l2 + 46 (aq) + 72 (aq) 2 + 56 (aq) 3 + 21 (aq) 4 + 3 (aqf^j / 

3 (6 + 12 (aq) + 9 (aq) 2 + 2 (aq) 3 ^ (2.13) 
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in which g (aq) is the correction to the classical result of Obreimoff (Obreimoff et 
al 1930) for peeling off a rigid substrate. 

Perturbation Analysis: Matching the coefficients for e l , i = 2,4 in the left 
and the right hand side of equation 2.11 results in the following set of equations: 

e 2 : Pix = U 1YY + U lzz , Piy = 0, Piz = 

e 4 : Pix = Uixx + ^2yy + Uizzi Piy = V\yy + Vizzi Piz = Wi YY + W\ zz 
U 1X + V 1Y + W 1Z = 

(2.14) 

we assume also that the excess quantities vary sinusoidally along the Y axis: so that 
Ti = fiSm(KY); T = U, V, W, P; i = 1, 2; K = 2tt/ (X/h) is the dimensionless 
wave number of the perturbed waves. Equations 2.14 are solved using the following 
boundary conditions derived from 2.9, 

(a) at Z = Ui = V x = W x = U 2 = 

(b) at Z = \ U lz = (U 2Z + W lx ) = (V lz + W 1Y )=0 (2.15) 

(c) at Z = 1,X<0 Pi=S^ixxxx, at < X < aq = ^ lxX xx 

where = W\ (X, Z = 1) is the vertical displacement of the flexible plate at Z = 1. 
Assumption of sinusoidal dependence on Y is a simplification which is apparent from 
the video micrographs of the contact line as in figure lb which show that the film 
and the plate remains in contact whole through the area of the finger implying 
that the deformation of the film is not perfectly sinusoidal, as it would mean a line 
contact between the plate and the film. However, here I assume sinusoidal variation 
to keep the calculations simple. Since, the plate does not bend in the direction 
of Y, and ^ixxxx both remain uniform along this axis; consequently in b.c. 
2.15c, bending stress on the plate is equated to Pi- Equation 2.14 suggests that 
Pi remains independent of Y and Z although U\ varies along Y, hence, the only 
solution for Pi that can satisfy equation 2.14 is Pi = 0. Other components of the 
excess displacements are obtained as: 




+ e- K ) \ 



~ K + 2K (e K - e- K ) , x 

Here, only the lowest order terms in e are computed since the higher order terms 
enhance accuracy insignificantly Since for all our experiments, e < 0.3, the above 
solutions imply that the excess deformations in the film occur under very small 
excess pressure which is of the order e 4 < 0.01. This excess pressure, however 
small, varies along Y implying that it should depend upon the distance between 
the plate and the film. Nevertheless, the excess traction which results from the 
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distance dependent forces (Shenoy et al 2001, Sarkar et al 2004) apply only in the 
immediate vicinity (< 0.1/zm) of the contact between the film and the plate as the 
gap between the two increases rather sharply (could be observed in AFM images of 
the permanent patterns of surface undulations). Hence, it does not contribute any 
significantly to the overall energetics. 

While equation 2.16 elaborates the variation of excess deformations along Y and 
Z, their dependence on X is incorporated through the coefficient C (X) which is 
obtained by solving equations 2.15c using the following boundary conditions 



(i. ii) *i|jr=o = c ^ (K) (iii) *i x \ x =o- = ^ix\ x =o+ 

(iv) *ixx\ x=0 - = *ixxl x=0 + ( v > vi ) *ilx=-« = *ixlx=- £ = 

(Vii, Viii) = VlXx\ X =a q = 



Where, £ = Aq is the dimcnsionless amplitude of the waves and Co$ (K) is the ex- 
cess stretching of the film at x = 0: Co is a constant and $ (K) — W\/C (X)\ x = 
(e~ K - e 3K + AKe K ) /K (l + e 2K (1 + 2K)). (a) B.C. i, ii, iii and iv occur because 
the excess displacement and slope of the plate are continuous at X = 0; (b) at 
X = — £ the displacement of the plate, its slope should vanish, which result in the 
b.c. v and vi; (c) similarly, at X = aq excess displacement and curvature of plate is 
zero, which results in b.c. vii and viii. Incorporating these boundary conditions into 
the solutions of equations 2.15c yield the following expression for the the excess 
displacement of the plate: 



^ = aq^iT+L) (" + 6ag$ + 2 {aq)2 ) ™ 3 " 3 + ^ ™ 2 

-3 (V 2 - 2 (aqf^j (X/£) + aq (3£ + 4aqfj at X < 

*i = C nK){l + ^^^ + 3aq) { X/aqf- 

3 (2£ + 3aq) {X/aqf - 3 ^ 2 - 2 {aqf^ (x/ (a?) 2 ))) at < X < aq (2.17) 

Excess Energy: Total energy of the system consists of the elastic energy of 
the film, bending energy of the plate and the intcrfacial energy: 



n = n e + n fe + n, 

r 2iv/k r h 

Uv z +w y ) + (uy + v x ) + (u z + w x ) J dzdydx 



oo JO JO 
Yy p a /*2Tr/k 



+— / (V>**) 2 dydx + W A (2na/k + A finger ) (2.18) 

* J -oo JO 

where Afi nger is the interfacial area of contact at —a < x < 0. From equation 2.18 
the excess energy of the system is obtained as H excess = U — U which is written 
in a dimensionless form using /j,/q 3 as the characteristic energy and by substituting 
for variables T = T + e 2 Ti + e 4 T 2 , where T = U, V, W and P. 
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The expression for excess elastic energy n e in the film is obtained as, 

n e = i y ? y o ^ (e 6 (V lz + W lY + e 2 (V 2z + W 2Y )) + 

e 8 (C/ 2y + V lx + e 2 V 2x ) 2 + e 8 {(U 2Z + W lx ) + e 2 W 2x f + 

2e 4 (J7oz + e 2 W 0x ) (U 2z + W lx + e 2 W 2x )) dZdYdX (2.19) 

where I estimate the excess energy within a distance — £ < X < 0. Considering only 
the leading order terms (e 4 and e 6 ) I can simplify the expression in 2.19 as, 

6 r r2n/K rl 

n e = - / / / (y lz + T^ ly ) 2 dZdrdx (2.20) 

4 «/-oo JO JO 

in which I substitute the expressions for V\ and Wi from 2.16 to obtain 

n e - e 6 C 2 f 2 ^,aq,K)/4 = e 6 C 2 f (K)f 2 ^,aq)/A (2.21) 
Similarly, the dimcnsionlcss excess bending energy lib of the plate: 

n b = (3ir/K) f q he 2 H> xx*ixx + e 4 (Hf lxx + e 2 ^ 2xx ) 2 ) dX (2.22) 

Considering only the leading order terms: e 2 and e 4 , equation 2.22 simplifies to 

n b = — / (2e 2 V Qxx y lxx +e^ 1 2 xx )dX (2.23) 

J — xi 

Substituting the expressions for an d ^1 from 2.12 and 2.17, lib is obtained as, 
n b = e 2 C Af 3 (£, aq, K) + e 4 C 2 fi (£, aq, K) (2.24) 
The interfacial energy is estimated as, 

n ; = 2Ge / | sin (KY) dY = — i- (2.25) 
Jo 2 K 

Substituting the expression for G from equation 2.13 in equation 2.25 and com- 
bining all the three energies, yields the total excess energy as: 

n = (e 6 / 2 & aq, K) + c 4 /i (£, 09, if)) C 2 + e 2 / 3 (£, aq, K) AC + ^^^ 2.26) 

The expressions for /i (£, ag, if), / 2 (£, aq, K) and / 3 (£, aq, K) are obtained using 
Mathematica and are not being presented here since they could not be written in 
a compact form. 



3. Results and Discussion 

The expression for excess energy in equation 2.26 accounts for the combined effects 
of three sets of parameters: e, the confinement parameter, £ and K, the characteris- 
tic length scale of the perturbations and aq and A, the length scale of the geometry 



Article submitted to Royal Society 



8 



A. Ghatak 



of the experiment. In what follows, we look for the solutions of these different sets of 
parameters which results in negative excess energy associated with the instability. 

While it is evident from equation 2.26 that £ and K are nonlinear ly coupled 
quantities, the physics of the problem is better understood if we study their effects 
separately. We do that following the observation that the excess displacements in 
the film but not that of the plate are functions of Y , which allows us to assume that 
the dimcnsionless wave number K of the perturbations is determined solely by the 
minima of the excess elastic energy II e and not the other components of the total 
excess energy. Although this assumption is not exactly correct as the displacement 
^ of the plate is also a function of K, experimental observation that the amplitude 
remains nearly independent of the wavelength (Ghatak et al 2003) suggests that 
the above assumption should not insert much inaccuracy into the calculation. 




1.4 1.6 1.8 ' 2.0 2.2 2.4 0.0 0.1 0.2 0.3 0.4 0.5 

K £ 

Figure 2. (a) Dimensionless elastic energy of the film is plotted against the dimensionless 
wave number K of the surface undulations. The energy of the film attains a minima at 
K — 1.91 implying that the wavelength varies with thickness of the film as A = 3.3h. 
(b)Amplitude £ are plotted against the limiting values of e from equation 3.2 for differ- 
ent values of the dimensionless length aq. Curves 1-4 represent aq — 5, 15, 25 and 55 
respectively. e u = 0.48 is the upper limit for e beyond which no real solution for Co exists. 



The plot of excess elastic energy II e vs. the wave number K in figure 2a then 
shows that II e attains a minima when K — 1.91. The wavelength of perturbations 
thus scales with the thickness of the film as A = 3.3/i which corroborates with 
the general observation in a wide range of experimental geometries (Ghatak et al 
2000, 2003, Monch et al 2001) that A remains independent of all the material and 
the geometric properties of the system except h. Furthermore, the proportionality 
constant matches well with that observed in experiments (A = 3 — 4/i)with rigid 
and flexible contacting plates. 

Although the minima of II e occurs at K = 1.91 irrespective of C and e, for 
both these parameters real values are desired. In fact, equation 2.26, being quadratic 
w.r.t. Co, suggests that in the limit II = the real solutions for Co exist only when, 

(e 2 f 3 (£, aq, K)f - 4 (e 6 / 2 (£, aq, K) + e\h (£, aq, K)) (27£/K) g (aq) / (aq) 4 > 0(3.1) 
resulting in the following inequality for e 2 : 

2 < h 2 (£> K ) K h (£, aq, K) 

" ~ Ah(Z,aq,K)27Zg(aq) f 2 ^,aq,K) l ' J 
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Equation 3.2 sets an upper bound for e as evident from figure 2b where £ and e 
which satisfy equation 3.2 are plotted for K = 1.91 and for different aq. When e 
is smaller than this upper critical limit e u , two different solutions for £ exist, the 
stability of which depends upon whether II attains a minima at these solutions. 
Hypothesizing that II minimizes when dYl/dC = 0, we obtain an expression for 
Co which when substituted in 2.26, yields 



n = - 



27A 2 £ 
H 7^9 (aq) 



(3.3) 



4 /i^a^ + ^iWO ' K 
In figure 3a we plot II from equation 3.3 w.r.t. £, for e = 0.48 — 0.05. For all these 



(a) 1.5 
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Figure 3. (a) Dimensionless excess energy is plotted against the dimensionless amplitude 
of the waves for different values of the confinement parameter e. The curves are obtained 
using representative values for the dimensionless parameters: aq — 25, and K = 1.91. (b) 
Bifurcation diagram showing variation of £ w.r.t. e for different values of aq and K — 1.91. 
The dotted and the solid lines represent respectively & vs. e and £ m i„ vs. e. No solution 
for £ exists beyond e c . 



cases, II exhibits a non-monotonic character: with increase in £, it first increases 
till it reaches a maxima after which it decreases to attain a minima, there after 
it increases again. Stability of these systems is determined by the minima of the 
excess energy which if positive signifies stable base solution for the contact line 
and unstable solution if negative. For example, for e = 0.48 II remains positive all- 
through implying that the undulation of a straight contact line would increase the 
total energy of the system so that a straight contact line would remain stable. On the 
other hand, e = 0.365, presents a limiting case for which minima of the excess energy 
n m j n attains zero and for e = 0.25, 0.15, 0.1 and 0.05, II m i n becomes negative, 
i.e. when the film is more than critically confined i.e. e < e c = 0.365, the contact 
line can become unstable if sufficiently perturbed. The critical value e c = 0.365 
thus obtained for the above set of data corroborates well with 0.3 obtained in 
experiments of figure 1 . Furthermore, a finite energy barrier at £ < & suggests that 
a straight contact line is not unstable for perturbations of all magnitude. Because, 
for perturbation with amplitude £ < & the excess energy remains positive so that 
these perturbations decay to zero. This result too corroborates with experiment that 
with increase in confinement of the film, the amplitude of the undulations never 



Article submitted to Royal Society 



10 



A. Ghatak 



increases from exactly zero, but from a finite value. When £ > II decreases 
to become negative, so that these perturbations can grow till £ reaches £ m , n , at 
which II attains the minima II m i n ; £ m i„ is then the predicted amplitude of the 
undulations of the contact line. 

Figure 3b depicts the bifurcation diagram where & (dashed line) and £ min (solid 
line) are plotted with respect to e for variety of aq. The dashed lines signify that 
the perturbations whose amplitude £ < £j decay to zero, whereas solid lines mean 
those with £ > & grow to £ TO j„. The amplitude £ mi „ increases with increase in 
the confinement of the film similar to that observed in experiments, although, the 
values predicted are somewhat ( 2-3 times) larger than what is observed. This 
discrepancy could be due to the underestimation of the excess elastic energy of the 
film. In figure 4 the combined effects of £ and K are probed by plotting II w.r.t. 




Figure 4. Excess energy II for e = 0.2 and aq = 25 is plotted w.r.t amplitude £ and wave 

number K. 

£ and K. Here again II m i n becomes negative for confinement parameter e below a 
critical value e c = 0.365. However, the wave number K at which the minima occurs 
does not remain constant, it decreases from 2.12 to 0.5 while e varies from 0.365 
to 0.05. Although this prediction is somewhat different from experiments, in which 
e varies between 0.3 to 0.1 at which K is observed to be 1.57 ± 0.1, some recent 
observations (personal communication with Prof. A. Sharma) with very thin elastic 
films (<~ 0.5/um) indeed indicate that K can decrease to 1.0 as e decreases to 0.07. 
More experiments are clearly necessary to characterize quantitatively the effect of 
the coupling of the two length scales with highly confined elastic films. 

4. Summary 

The analysis shows that confinement of an incompressible elastic film leads to favor- 
able energetics for perturbations to grow so that the film can not deform uniformly 
everywhere when subjected to the tensile stresses at the interface. Furthermore, 
the nature of the adhesion stress is not important, even the spatial variation of 
the surface forces play rather an insignificant role. While the theory captures the 
essential physics of the problem, slight overestimation of the amplitude possibly 
results from the assumption of sinusoidal variations along the Y axis which is not 
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perfectly correct. These issues can possibly be resolved by 3D simulation of the 
force field near the contact line. Nevertheless the results presented in this paper 
should be important for many other systems in confined geometries. 
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